Hierarchy of fractional Chern insulators and competing compressible states 



A.M. Lauchli/'^ Z. Liu,^ E.J. Bergholtz,"^' ^ and R. Moessner^ 

^Institut fur Theoretische Physik, Universitdt Innsbruck, A-6020 Innsbruck, Austria 
^ Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzer Strafie 38, D-01187 Dresden, Germany 
^Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China 
^Dahlem Center for Complex Quantum Systems and Institut fur Theoretische Physik, 
Freie Universitdt Berlin, Arnimallee 14, 14195 Berlin, Germany 
(Dated: July 26, 2012) 

We study the phase diagram of interacting electrons in a dispersionless Chern band as a function of their 
filling. We find hierarchy multiplets of incompressible states at fillings = 1/3, 2/5, 3/7, 4/9, 5/9, 4/7, 3/5 
as well as f = 1/5, 2/7. These are accounted for by an analogy to Haldane pseudopotentials extracted from 
an analysis of the two-particle problem. Important distinctions to standard fractional quantum Hall physics are 
striking: absent particle-hole symmetry in a single band, an interaction-induced single-hoX^ dispersion appears, 
which perturbs and eventually destabilizes incompressible states as v increases. For this reason the nature of the 
state at = 2/3 is hard to pin down, while zy = 5/7,4/5do not seem to be incompressible in our system. 

PACS numbers: 73.43.-f, 71.10.Fd, 73.43.Nq, 



Introduction. — Following recent proposals of the existence 
of novel lattice generalizations of fractional quantum Hall 
(FQH) states, termed fractional Chern insulators (FCI), in 
(approximately) flat bands exhibiting non-zero Chern num- 
bers [1-3], there has been intense research activity in under- 
standing this phenomenon [4-26]. 

From the original observation of a FCI state at Chern band 
filling = 1/3 [3-5] (at which the original FQH state was 
also first observed), a number of questions immediately arise. 
Firstly, under what conditions can FCIs be observed? Sec- 
ondly, what are the states which compete with the FCI states? 
Thirdly, what are the differences between FCI physics in 
Chern bands compared to the familiar setting of Landau levels 
in the continuum appropriate for describing the FQH state in 
conventional semiconductor heteroj unctions, e.g. arising due 
to the non-uniform Berry curvature in reciprocal space? 

This publication aims to contribute to all of these questions. 
We start by demonstrating that a nearest neighbor interaction 
leads to FCI state also dXv = 2/5 and u = 1/5, although the 
latter fraction has an FCI phase that is substantially less robust 
than the = 1/3 FCI. We account for this with a heuristic de- 
rived from considering the two-particle problem in the lattice 
model. Further, we find evidence of several additional FCI 
states akin to the hierarchy FQH states familiar from conven- 
tional QH physics and its composite Fermion [27] hierarchy 
[28, 29] picture. 

Turning to the qualitative distinctions from conventional 
QH physics, we find that in the absence of particle-hole 
symmetry in a single band, an effective crystal-momentum- 
dependent potential appears, which shows up in the properties 
of the single-hole dispersion, and leads to a modulation of 
the occupation numbers n(k) in the many-body ground state 
which presage the breakdown of the QH effect. This is re- 
flected by a strong distinction between the many-body states 
atu = 1/3 and u = 2/3 even in the spin-polarised setting that 
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Figure 1. (Color online) (a) Illustration of the checkerboard lattice 
model with the relevant hopping amplitudes [2, 3]. The interaction 
term is a nearest neighbor density-density interaction (interaction be- 
tween pairs of white and grey sites), (b) Berry curvature in the first 
Brillouin zone for the indicated set of parameters. Dark (light) inten- 
sity denotes small (large) Berry curvature. 



we study. 

Our studies use extensive large-scale exact diagonalisations 
of the many-body Hamiltonian of the lattice systems coupled 
with analytical considerations. On a more technical level, we 
address questions of the interplay of topological order and 
the concomitant finite- size lattice quasi-ground- state degen- 
eracies, as well as the finite-size scaling of the gaps. 

Setup. — In order to be specific we focus on the two- 
band checkerboard lattice model introduced in Refs. [2, 3] 
and sketched in Fig. 1(a). Here tie^*"^ is a nearest-neighbor 
hopping with an orientation dependent complex phase, and t2 
denotes the next-nearest-neighbor hopping amplitude. After 
Fourier transform, the single-electron (kinetic) Hamiltonian 
reads 

n=J2{cl^,cls)Hk){c^A,c^BV, (1) 

kGBZ 

where /i(k) is given as [36]. In the following we con- 
sider the case ofti = 1,(/) = 7r/4 and varying ^2. For 
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Figure 2. (Color online) (a) Energy of a single hole in the fully filled 
lower band ^^(K) := E{Np = iV, - 1, K) - E{Np = A^^, 0). (b) 
Non-zero eigenvalues E{Np = 2, K) of the two particle problem 
on a representative path through the Brillouin zone for t2 = 0.5. In 
general there are four non-zero energies except at K = (0,0) where 
there are only two finite energies. Circles are data for (b) 6 x 6 (empty 
circles) and 8x8 (filled circles) unit cells [(a) 10 x 10 unit cells], 
and the lines are guides to the eye. 

^2 7^ 0, oo one obtains two separated bands with Chern num- 
ber ±1. The Berry curvature in each band can be expressed 
in terms of the single-particle states in the band, via 
the Berry connection ^^(k) = — ^(^k l^^j l^k ) -^ij ~ 
9/c.A^(k) — dkjAfCk). The two bands have topologically 
quantized Chern numbers, C = ^ F^{k)d?k = ±1, 
while the Berry curvature varies in the Brillouin zone and de- 
pends on the microscopic parameters (t2 in the present case). 
In Fig. 1(b) we display the Berry curvature for the specific 
value of ^2 = 0.5 [note that the Berry curvature vanishes at 
(0,0) and (7r,7r)]. 

Interactions. — We consider nearest neighbor repulsion, 
^int = [cf. caption of Fig. 1(a)], throughout 

this publication. If the energy gap between the two bands 
is large compared to the interaction strength it is possible to 
project the interactions into the partially filled band, leading 
to a projected Hamiltonian of the form 

H= ^ Vkik2k3k4Cki42^k3Ck4, (2) 

kik2k3k4 

for two-body interactions. Here, cq annihilates an electron 
with wavefunction ^Q(r) of the single particle state of the 
partially filled band with crystal momentum Q. In the follow- 
ing we take the flat band limit, i.e. drop the residual energy 
dispersion in the band, as we wish to highlight the new phe- 
nomena arising due to the non-constant Berry curvature, as 
opposed to remnant energy dispersion. 

Particle-hole asymmetry and two -particle analysis. — It is 
enlightening to consider the two simplest interacting cases, 
namely the single-hole and the two-particle problem. Per- 
forming a particle-hole transformation, Ck c\. within one 
of the bands, the projected Hamiltonian (2) transforms to 

J2 ^k.k.k3k4«Ck3Ck4 + E^''(l^)4ck(3) 

kik2k3k4 k 

which includes an effective single-hole energy £^^(k) = 
4 Xlm ^mkmk (adopting a convention under which Vkimn 



is antisymmetric under exchange of its indices). This term 
amounts to a trivial overall energy shift in a Landau level. 
By contrast, it introduces an effective dispersion even for an 
entirely flat Chern band! We display E^ili) in Fig. 2(a) for 
t2 = 0.5. As we will demonstrate below, this particle-hole 
asymmetry generally leads to a deformation of n(k) and in 
fact dominates the physics near v = 1, i.e. in the dilute-hole 
limit. Note that for our checkerboard example, the energy 
spectrum is particle-hole invariant under u ^ 2 — v. 

Next, we consider the two particle problem: In Fig. 2(b) 
we display the non-zero eigenvalues of the two particle prob- 
lem along a path in the Brillouin zone for the nearest neigh- 
bor interaction for ^2 = 0.5. The spectrum depends on 
the total momentum, K = ki + k2, underscoring the lack 
of translation invariance in reciprocal space. We observe 
two dominant eigenvalues of mean value ^ 0.4F which re- 
main non zero throughout the Brillouin zone, while there is 
also second set of two eigenvalues about 10 times smaller 
which vanish as K (0,0). All other eigenvalues are 
strictly zero. Quite generally, for each K, the number 
of non-zero energy levels is bounded above by the num- 
ber of finite energy levels of the interaction alone before 
band projection. This follows from the fact that each (un- 
projected) interaction term imposes one (linear) constraint 
which needs to be satisfied for the two particle wavefunction 

IV^k) = Eq^q [</>K+q(ri)(/>-q(r2) - 0K+q (^2 ) ^-q (ri )] tO 

have zero energy: (V^k| V^(^i, ^2)|V^k) = 0. If these con- 
straints are linearly dependent, there can be additional zero- 
energy states, as is the case at K = (0, 0) in the present exam- 
ple. Here, the upper bound of 4 finite eigenvalues is saturated 
everywhere else in the Brillouin zone. Including also a next- 
nearest neighbor repulsion leads to 8 finite eigenvalues, and 
e.g. on a three- sublattice kagome lattice model [1] one finds 6 
(12) levels including the (up to next-)nearest neighbor terms. 

These observations motivate a suggestive analogy between 
the finite energy levels of the projected two-particle problem 
and Haldane's pseudopotenials well known in the continuum 
Landau levels [28]. We thus tentatively label the two non-zero 
pairs of eigenvalues as "pseudopotentials" Vi and V3. 

Hierarchy states and the phase diagram. — To investigate 
the full interacting many-body problem we have performed 
extensive exact diagonalization studies on a large number of 
finite samples with rectangular shapes (as in previous numer- 
ical studies of FCI) as well as tilted samples, where the span- 
ning vectors of the samples need not be aligned with the lattice 
axes. The latter choice of samples allows to study a consider- 
ably larger number of clusters with aspect ratio close to one. 
The list of considered clusters can be found in the supplemen- 
tary material [37]. 

The pseudopotential analogy introduced above suggests a 
hierarchy of incompressible states dX v = p/q > 1/3, q 
odd, due to the relatively large energy scale Vi. In particu- 
lar, we have looked for the m = 1 composite fermion branch 
at ly = p/ {2mp + 1) [27] with rather clear evidence of FCI 
states at ly = 2/5,3/7, 4/9 as well as their 1 — u relatives 
aiu = 5/9, 4/7, 3/5, beyond the u = 1/3 state already dis- 
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Figure 3. (Color online) (a)-(c) Numerical evidence for a Fractional Chern insulator state at u = 2/5. For details see main text, (d)-(f) 
Numerical evidence for a Fractional Chern insulator state at = 1/5. Here we focus on t2 = 1/2 which is close to the value where the 
variance of the Berry curvature is minimal. 



cussed in the literature [3-5]. Data of the low-energy spec- 
trum of at least two samples for each fraction are displayed in 
the supplementary material [37]. Here we provide strong ev- 
idence for a stable incompressible 1/ = 2/b state as shown in 
Fig. 3(a)-(c). In panel Fig. 3(a) we present the spectral flow of 
the 5 ground state manifold levels upon flux insertion for two 
different samples. The ground state manifold (filled symbols) 
does not mix with excited states (empty symbols) and returns 
to the initial spectrum upon insertion of five flux quanta. This 
provides strong evidence for a state with quantized Hall con- 
ductance. In Fig. 3(b) we display the finite size scaling of the 
energy spread of the five states in the ground state manifold 
for different samples. We plot the energy splitting as a func- 
tion of the inverse topological diameter 1/W (sls defined in the 
supplementary material [37]), which can become quite small 
(as small as 1/25) for our tilted samples. Note that the topo- 
logical extent is not in general equivalent to the geometrical 
extent. In particular the samples with W = 10, 20, 25 have 
a geometrical aspect ratio equal to one. We also stress that 
the five ground states are systematically found in the sectors 
predicted by the counting rule developed in Refs. [5, 9], and 
which we generalized for tilted samples [37]. In the case of 
tilted samples the ground state sectors show a great variabil- 
ity depending on the cluster geometry, which is an argument 
against the formation of a charge density wave state, which 
should exhibit a unique set of degenerate momenta dictated 
by the spatial symmetry breaking of the charge density wave. 



Finally, in Fig. 3(c) we show the finite size scaling of the gap 
from the ground state to the first excited state as a function 
of inverse systems size 1/Ns. Our data shows non-negligible 
finite size effects, but convincingly points towards a finite ex- 
citation gap of the order of V/ 30 in the thermodynamic limit, 
in contrast to recent claims denying a stable u = 2/b state for 
the checkerboard lattice [16, 19]. 

Next we turn to fractions below u = 1/3. According to 
the pseudopotential analogy our nearest neighbor interactions 
generate a finite V3, albeit roughly a factor ten smaller than 
Vi. We would thus expect to find a stable u = 1/5 state 
with associated energy scales significantly smaller than those 
of the = 1/3 state. Indeed we find strong evidence for a 
v = 1/5 FCI state as shown in Fig. 3(d)-(f). In particular the 
finite size effects of the ground state manifold splitting shown 
in Fig. 3(e) behave qualitatively similar to the = 2/5 case 
in Fig. 3(b), and the energy gap to excited states extrapolates 
very nicely to a value of the order of V/200. We furthermore 
find some mild evidence for a = 2/7 state [37], which is 
stabilized by V3. 

Effect of the single hole energy. — When numerically ex- 
ploring the phase diagram we discovered that several particle 
hole conjugate states of stable low-density FCI fractions do 
not seem to be realized, for example we find it difficult to sys- 
tematically observe the required ground state degeneracy for 
u — 2/3, while for v = 5/7 and u = 4/5 the required de- 
generacy is absent. An enlightening way to understand this 
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Figure 4. (Color online) left panel: Momentum space orbital oc- 
cupation n(k) = (c^Ck) plotted as a function of the single-hole 
energy Eh(k.). As the filling u is increased, n(k) displays a more 
pronounced (close to linear) correlation with ^^(k). Dashed lines 
denote fractions which are most likely not FCI states. Right panel: 
Approximate slope of the correlation between n(k) and Eh (k) plot- 
ted as a function of the filling v. The rapid, nonlinear increase of 
the slope for large fillings is a likely cause for the breakdown of FCI 
states at large filling. 



finding is to monitor the filling dependence of the momentum 
space occupation number n(k) = (cj^Ck) when plotted as a 
function of the single hole energy Eh(k) introduced before, 
see Fig. 4(a), where the occupancy n(k) in the incompress- 
ible many-body state clearly tracks £^/i(k). The role of the 
interaction-induced effective single particle dispersion is thus 
to significantly distort the momentum space occupancy as the 
filling increases, and ultimately leads to Fermi- surface like 
compressible states which are governed by the effective sin- 
gle particle dispersion. This behavior is illustrated in Fig. 4(b) 
with the drastic increase of the correlation between n(k) and 
EhCk) as ly increases, and the Fermi-surface like n(k) in the 
z/ = 4/5 state depicted in Fig. 4(a). 

It is instructive to compare our results from what is known 
from conventional FQH physics. For intermediate filling frac- 
tions in the lowest Landau level, both experiment [30] and the- 
ory [31] suggest that the stability of the FQH states towards 
disorder, temperature etc., is essentially determined by the de- 
nominator, q, of the filling fraction (for odd q) implying a self- 
similar structure in the clean zero-temperature limit. Correc- 
tions to the self-similar structure in "real" systems are minor 
and include effects due to an instability towards Wigner crys- 
tallization at low filling fractions and particle-hole symmetry 
breaking effects (LL mixing etc). The main features are qual- 
itatively similar to what we find in the lattice system. How- 
ever, in the Chem insulator case the corrections are more sub- 
stantial. In particular, there is nothing like a "clean" limit as 
the Berry curvature effects necessarily destroys weak states, 
and the compressible competing states discussed above pre- 
vail. Moreover, due to the explicit particle-hole asymmetry 
within the band we find that FCI states are absent for z^ > 2/3. 



Discussion. — The question to what extent fractional 
Chern insulators are 'nothing but' FQH states on the lattice 
is fundamental to this rapidly developing field: it combines 
the questions of where to find a FCI with how to diagnose it 
and how to theoretically describe it, before being able to ad- 
dress differences to the 'weak-field' continuum Landau level 
physics. 

The FCI states we thus found, as well as their respective 
stabilities, are naturally organised with the aid of the hierar- 
chy concept developed for the FQH. To obtain this informa- 
tion, we have developed new diagnostics, including an adap- 
tation of Haldane's pseudopotentials, strongly suggesting a 
common field-theoretic description for the two phenomena. 
At the same time, a particularly noteworthy negative result is 
our observation, at odds with some previous treatments, that 
FCIs do not generically come with a constant occupation num- 
ber of the flat band orbitals in reciprocal space. 

Indeed, this phenomenon arises even for perfectly flat (dis- 
persionless) bands, as a result of a non-uniform Berry curva- 
ture, which appears as a potent driving force for the destruc- 
tion of FCI states by acting as a /c-dependent chemical poten- 
tial of a strength tuned by the electron-density of the flat band, 
removing particle-hole symmetry in a lattice- specific manner. 
Our study provides a a first glimpse of the competing states. 

Clearly, much further work on each of these points is 
warranted, especially in the light of promising experimental 
prospects to observe such states in optical lattices with artifi- 
cial magnetic fields [20, 33-35]. 

Note added. — While completing the present work we be- 
came aware of two preprints which also provide evidence for 
FCI states beyond the Laughlin states [16, 19]. 

We acknowledge useful discussions with Masud Haque. 
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"Scientific computing" at the University of Innsbruck - sup- 
ported by the BMWF - and on the PKS-AIMS and vip cluster 
at the MPG RZ Garching. 
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Figure 5. (Color online) Geometry of a A^s = 20 site cluster with 
the spanning vectors Ti and T2, its momentum space point and an 
illustration of the topological extent. The shaded grey momentum 
points correspond to the expected ground state momenta for = 1 /5 
state derived by the counting rule for a 10 x 2 lattice. 





Spanning vectors 


Topological extent 




T1/T2 


W xU 


18 


(3,3)/(-3,3) 


6x3 


20 


(4,2)/(-2,4) 


10 X 2 


25 


(5,0)7(0,5) 


5x5 


25 


(4,3)/(-3,4) 


25 X 1 


28 


(4,0)7(0,7) 


4x7 


30 


(5,0)7(0,6) 


5x6 


35 


(5,0)7(0,7) 


5x7 


36 


(6,0)7(0,6) 


6x6 


40 


(6,2)7(-2,6) 


20 X 2 


45 


(6,3)7(-3,6) 


15 X 3 


50 


(5,5)7(-5,5) 


10 X 5 


50 


(7,1)7(-1,7) 


50 X 1 



Table I. List of samples studied in the present work. Apart from 
the rectangular samples with 28, 30 and 35 sites, all clusters have an 
geometric aspect ratio of 1 and at least a C4 point group symmetry 
in addition to the translation group. 

SUPPLEMENTARY MATERIAL 
Tilted Square samples 

We use rectangular and tilted samples in our exact diago- 
nalization study. The samples are defined by their spanning 
vectors Ti and T2, which define the toroidal periodic bound- 
ary conditions. A list of clusters studied in this work is shown 
in Tab. L The tilted clusters used here all have perpendicular 
spanning vectors Ti and T2 of equal length, but these require- 
ments could also be relaxed in future studies for an even larger 
sample variability. An interesting aspect of the tilted samples 
is their so called topological extent. While they are symmetric 
samples of aspect ratio one in real space [cf. Fig. 5(a)], they 
actually have a different aspect ratio in momentum space. We 
define the topological extent W as the number of steps one 
has to go in momentum space on a straight path along nearest 
neighbor orbitals in momentum space until one returns to the 
starting point. The complementary extent U is then simply 
defined as [/ = Ns/W. The topological extent of all clus- 



ters is listed in Tab. L The extents W and U are important 
for two reasons: i) it seems that the parameter W controls 
to a large extent the splitting of the ground state manifold as 
demonstrated for z/ = 1/5 and u = 2/5 in Fig. 3(b), (e) of 
the main paper. This then allows to probe for the existence of 
a degenerate ground state manifold on much smaller clusters 
compared to the usually employed rectangular clusters, where 
the topological extent is identical to the real space extent, and 
thus values for W of the order of ten or more will be difficult 
to reach with exact diagonalization. The second reason ii) is 
that the counting rule for ground state quantum number pre- 
diction [5, 9] can then be applied to the topological W x U 
lattice, and then be embedded back into the original momen- 
tum space to obtain the expected ground state momenta at a 
certain filling. This is illustrated by the shaded circles in Fig. 5 
(b) and (c) for a specific filling fraction. 



Numerical support for hierarchy states 

In the following we present low energy spec- 
tra for various fraction where there is some ev- 
idence for FCI states. These are the fractions 
u = 1/5,2/7,1/3,2/5,3/7,4/9,5/9,4/7,3/5. We 
find the situation at z/ = 2 /3 to be unclear, while v = 5/7 
and ly = 4/5 do not show any evidence for FCI states. For 
the bold fractions we now show the low energy spectrum for 
two or more samples. For the FCI fractions the ground state 
exhibits the correct approximate degeneracy at the momenta 
predicted by the counting rule [5, 9] in some range of ^2 
values. 

Apart from the fraction 1/ = l/b and u = 2/b discussed 
in the main paper (as well sls u = 1/3 discussed in the liter- 
ature [3-5]) it is quite difficult to perform a systematic finite 
size scaling study, because of the sometimes pronounced sam- 
ple dependent results, as well as due to limited sample avail- 
ability for certain fractions. In all of the following plots the 
expected ground state sectors are denoted by the filled sym- 
bols in different colors, whereas the other sectors are denoted 
by hashed symbols. The range of t2 values where the cor- 
rect approximate ground state degeneracy can be found be- 
fore the first excited state is denoted by an orange box in 
each graph. As the data is shown it is in principle only evi- 
dence for the emergent many-body translation symmetry, sim- 
ilar to the FQH states on the torus, where any state at filling 
v = p/q is g'-fold degenerate. We however also checked the 
flux threading property of the ground state manifold - similar 
to Fig 3(a),(d) - for several of the fractions, with a positive 
results. One could also perform quasihole counting, but this 
requires care because of aliasing problems, for example on the 
Ns = 35 site sample the states at = 2/5 and u = 3/7 (and 
at = 3/5 and u = 4/7) just differ by one electron. 
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Figure 6. (Color online) FCI at z/ = 2/7 
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Figure 7. (Color online) FCI at = 2/5 



8 



< • (0,0) 

^ H (1,0) (4x) 
(2,0) (4x) 
A (1,1) (4x) 
< (2,1) (8x) 
V (2,2) (4x) 



v=3/5, Ng=25 (5,0)/(0,5) 



A 







i: 



A 
4 



V 




v=3/5, N =35 (5,0)/(0,7) 




• 


(0,0) (1x) 


• 


(1,0) (2x) 


• 


(2,0) (2x) 




(0,1) (2x) 




(1,1) (4x) 




(2,1) (4x) 




(0,2) (2x) 




(1,2) (4x) 




(2,2) (4x) 


A 


(0,3) (2x) 


A 


(1,3) (4x) 




(2,3) (4x) 



Figure 8. (Color online) FCI at = 3/5 
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Figure 9. (Color online) FCI at = 3/7 and = 4/7 
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Figure 10. (Color online) FCI at = 4/9 and = 5/9 
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Figure 11. (Color online) = 2/3. Unclear situation regarding a FCI state. 
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Figure 12. (Color online) u = 4/5. No evidence for a FCI state. 



